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Abstract 

We present expressions for shape variables of B decay distributions in several different mass 
schemes, to order a'^Po and Aqq-^/^I- Such observables are sensitive to the b quark mass and 
matrix elements in the heavy quark effective theory, and recent measurements allow precision 
determinations of some of these parameters. We perform a combined fit to recent experimen- 
tal results from CLEO, BABAR, and DELPHI, and discuss the theoretical uncertainties due to 
nonperturbative and perturbative effects. We discuss the possible discrepancy between the OPE 
prediction, recent BABAR results and the measured branching fraction to D and D* states. We 
find \Vcb\ = (40.8 ± 0.9) x 10-^ and mj'^ = 4.74 ± 0.10 GeV, where the errors are dominated by 
experimental uncertainties. 
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I. INTRODUCTION 



The study of flavor physics and CP violation is entering a phase when one is searching 
for small deviations from the standard model. Therefore it becomes important to revisit the 
theoretical predictions for inclusive decay rates and their uncertainties, which provide clean 
ways to determine fundamental standard model parameters and test the consistency of the 
theory. 

Experimental studies of inclusive semileptonic and rare B decays provide measurements 
of fundamental parameters of the standard model, such as the CKM elements \ Vcb\-, \Vub\-, and 
the bottom and charm quark masses. Inclusive and rare decays are also sensitive to possible 
new physics contributions, and the theoretical computations are model independent. The 
operator product expansion (OPE) shows that in the ^ Aqcd limit inclusive B decay 
rates are equal to the h quark decay rates [1, 2], and the corrections are suppressed by 
powers of and Aqcd/^^?)- High-precision comparison of theory and experiment requires a 
precise determination of the heavy quark masses, as well as the matrix elements Ai 2, which 
parameterize the nonperturbative corrections to inclusive observables at (9(Aqcd/^6)^- At 
order (Aqcd/''^6)^, six new matrix elements occur, usually denoted by pi,2 and ^,2,3,4- There 
are two constrains on these six matrix elements, which reduces the number of parameters 
that affect B decays at order (Aqcd/''^&)^ to four. 

The accuracy of the OPE predictions depends primarily on the error of the quark masses, 
and to a lesser extent on the matrix elements of these higher dimensional operators. It was 
proposed that these quantities can be determined by studying shapes of B decay spec- 
tra [3-6]. Such studies have been recently carried out by the CLEO, BABAR and DELPHI 
collaborations [7-13]. A potential source of uncertainty in the OPE predictions is the size 
of possible violations of quark-hadron duality [14]. Studying the shapes of inclusive B de- 
cay distributions may be the best way to constrain these effects experimentally, since it 
should influence the relationship between shape variables of different spectra. Thus, test- 
ing our understanding of these spectra is important to assess the reliability of the inclusive 
determination of \Vcb\^ and also of \Vub\- 

In this paper we present expressions for lepton and hadronic invariant mass moments 
for the inclusive decay B — > XJ,D, as well as photon energy moments in S — > decays. 
We give these results as a function of cuts on the lepton and photon energy, respectively. 
Most results in the literature have been given in terms of the pole mass, which introduces 
artificially large perturbative corrections in intermediate steps, making it difficult to estimate 
perturbative uncertainties. We present all results in four different mass schemes: the pole 
mass, the IS" mass, the PS mass, and the MS mass. We then carry out a combined fit 
to all currently available data and investigate in detail the uncertainties on the extracted 
parameters |Vc5| and nib. 

The results of this paper can be combined with independent determinations of the h and 
c quark masses from studies of QQ states. We have chosen not to discuss those constraints 
here, since there exist many detailed analyses in the literature [15]. Furthermore, the de- 
terminations of iTib and rric from QQ states have theoretical uncertainties which are totally 
different from the current extraction. Consistency between the extractions is therefore a 
powerful check on both determinations. 
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II. SHAPE VARIABLES 



We study three different distributions, the charged lepton spectrum [3, 4, 16, 17] and 
the hadronic invariant mass spectrum [5, 16, 18] in semileptonic B XJ,D decays, and the 
photon spectrum in S — > Xg'^ [6, 19-21]. Similar studies are also possible in S ^ X^f+f" 
and B XsVP decay [22], but at the present, such processes do not give competitive 
information. 

The B Xc(-i> decay rate is known to order q;^/?o [23] and Aq^j^/m^ [16], where /3o = H — 
2n//3 is the coefficient of the first term in the QCD /3-function, and the terms proportional 
to it often dominate at order a^. For the charged lepton spectrum we define shape variables 
which are moments of the lepton energy spectrum with a lepton energy cut. 



/ 



Kol^-C/O, -C/lj — — -rp , Kn[J^o) — -Tp , UJ 

^dEe / -—dEe 

Jeo dEi JEo dEi 

where dV/dEi is the charged lepton spectrum in the B rest frame. J?„ has dimension GeV"^, 
and is known to order a^/?o [17] and AQQ-^/ml [16]. Note that these definitions differ slightly 
from those in Ref. [4], and follow the CLEO [9, 10] notation. The DELPHI collaboration [12] 
measures the mean lepton energy and its variance (both without any energy cut), which are 
equal to -Ri(O) and -R2(0) — i?i(0)^, respectively. 

For the B — * Xciiy hadronic invariant mass spectrum we define the mean hadron invariant 
mass and its variance, both with lepton energy cuts Eq, 

Si{Eo) ^{ml-ml)_^, S,{Eo) = ({ml - {m\)f) _ ^ , (2) 



Ei>Eo 



Ee>Eo 



where m/j = {rriD + 3mD*)/4 is the spin averaged D meson mass. It is conventional to 
subtract rnjy in the definition of the first moment Si{Eq). Sn has dimension (GeV)^", and 
is known to order a'^/3o [18] and Aq^j^/mf [16]. For a given Eq, the maximal kinematically 



allowed hadronic invariant mass is = \Jm% — 2mBEQ. Once rn^^'^ — rfiD ^ Aqcd, the 
OPE is expected to describe the data. 

The above shape variables can be combined in numerous ways to obtain observables that 
may be more suitable for experimental studies because of reduced correlations. For example. 
Si and Rq can be combined to obtain predictions for 

, 2__2. _ Si{E^)R,{{),Eo) - Si{Ei)Ro{0,Ei) 

{mx m^)^^^^^^^^- R,(o,Eo)-Ro{0,Ei) ' ^"^^ 

that allows comparing regions of phase space that do not overlap [24] . 

For B Xgj, we define the mean photon energy and variance, with a photon energy cut 

Eq: 

TiiEo) = (E,) ^ UEo) = ({E, - {E,)f) _ ^ , (4) 



E-,>Ea' 



where dV/dE^ is the photon spectrum in the B rest frame. Again, T12 arc known to order 
alPo [19] and AQ(-;j-,/m^ [21]. In this case the OPE is expected to describe Ti(£'o) once 
itib/^ — Eo ^ Aqcd- Precisely how low £'0 has to be to trust the results can only be decided 
by studying the data as a function of £^0; one may expect that Eq — 2GeV available at 
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present is sufficient. Note tliat tlie perturbative corrections included are sensitive to tlie 
mc-dependence of tlie b — >■ ccs four-quark operator (O2) contribution. This is a particularly 
large effect in the O2 — O7 interference [19], but its relative influence on the moments of the 
spectrum is less severe than that on the total decay rate. The variance, T2, is very sensitive 
to any boost of the decaying B meson; this contribution enhances T2 by /9^/3 at leading 
order [19], where f3 is the boost {f3 ^ 0.064 if the B originates from T(4S') decay). This is 
absent if dr/d£^^ is reconstructed from a measurement of dV/dE^-^^. 



III. MASS SCHEMES 

The OPE results for the differential and total decay rates are given in terms of the b 
quark mass, m^, and the quark mass ratio, rric/mb. (Throughout this paper quark masses 
without other labels refer to the pole mass.) The pole mass can be related to the known 
meson masses via the l/mq expansion 

, T Ai + dMA2(mQ) Pi + dMP2 Ti + % + dui^ + %) , 

"^- = "^^+^ — + ^^^5 

where ruM (M — P, V) is the hadron mass, tuq is the heavy quark mass, and dp — 3 for 
pseudoscalar and dy — —I for vector mesons. The Aj's and p^'s are matrix elements of local 

dimension-5 and 6 operators in HQET, respectively, while the T^'s are matrix elements of 
time ordered products of operators with terms in the HQET Lagrangian, and arc defined 
in [16]^. The ellipses denote ^qcb/^q corrections, which can be neglected to the order we 
are working. Using Eq. (5), we can eliminate rric in favor of and the higher order matrix 
elements, 

m,-m. = m.-m,-A,(2^-^) + (p,-r,-r3)(^-^), (6) 

where Ttim = {mp + 3my)/4 denotes the spin averaged meson masses. 

Only three linear combinations of 7^_4 appear in the expressions for B meson decays (a 
fourth hnear combination would be required to describe B* decays). The reason is that the 
7i_4 terms originate from two sources: (i) the mass relations in Eqs. (5) and (6) which depend 
on Ti + T^ and T2 + T4; and (ii) corrections to the order Ag^p/m^ terms in the OPE, which 
amount to the replacement Ai Ai + (Tl + 372) / nib and A2 — >■ A2 + {% + ST^) / {Srub) . Since 
T1 + 3T2 = (71 + 7^) + 3(7^ +74) — (7^ + 37^), only three linear combinations are independent. 
Therefore, we may set 74 = 0, and the fit then projects on the linear combinations 

Ti-3r4, T2 + T4, r3 + 3T4. (7) 

The mass splittings between the vector and pseudoscalar mesons, 

AruM ^my-mp^ ^ "^"^^^ ^'^"^'^ - " + + . . . , (8) 

niQ 



^ These are related to the parameters p^,, p^^, p^^, p^Q, p% and p\ introduced in [25]. 
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constrain the numerical values of some of the HQET matrix elements. Here K{mc) — 
[as{mc) /as{'mh)]^^'^° ~ 1.2 is the scaling of the magnetic moment operator between mi, 
and rric. In terms of the measured B* — B and D* — D mass splittings, Ams and Am^, 

\ (m) - ^fe - ml Amp 
^ 2[mb - K{mc) rric] 

P2-(^+^) = 7 ^ ■ (10) 

These equations differ slightly from those in Ref. [16], and are consistent to order I/itlq. 
Since order Q;<i(AQCD/''77(5)^ corrections in the OPE have not been computed, whether we set 
K{mc) to its physical value, n{mc) — 1.2, or to unity is a higher order effect that cannot be 
consistently included at present. Using K{mc) — 1.2 or 1 in the fits gives effects which are 
negligible compared with other uncertainties in the calculation. 

It is well-known that the pole masses suffer from a renormalon ambiguity [26], which 
only cancels in physical observables against a similar ambiguity in the perturbative expan- 
sions [27]. Although any quark mass scheme can be used to relate physical observables to 
one another, the neglected higher order terms may be smaller if a renormalon-free scheme 
is used. When using pole masses it is important to always work to a consistent order in 
the perturbative expansion, since A can have large changes at each order in perturbation 
theory, even though the relations between measurable quantities such as the shape variables 
and the total semileptonic decay rate have much smaller changes. Since A depends strongly 
on the order of the calculation in perturbation theory, one can get a misleading impression 
about the convergence of the calculation, and its uncertainties. The advantage of using 
renormalon-free mass schemes is that the convergence may be manifest. 

Several mass definitions which do not suffer from this ambiguity have been proposed in 
the literature, and wc consider here the MS, IS*, and PS masses. (There is a renormalon 
ambiguity in the IS* and PS masses, but it is of relative order Aq^p/m^ and so is irrelevant 
for our considerations.) The MS mass is related to the pole mass through 

^^ = l-e^ii^ -6^1.562 ^#£/3o + .... (11) 

777.6 TT TT 

and Cp — 4/3 in QCD. The parameter e = 1 is a new expansion parameter, which for 
the MS mass is the same as the order in a^. While the MS mass is appropriate for high 
energy processes, such as Z or /i — > 66, it is less useful in processes where the typical 
momenta are below 7?76. The MS mass is defined in full QCD with dynamical 6 quarks and 
is appropriate for calculating the scale dependence above 7?76. However, it does not make 
sense to run the MS mass below 7775; this only introduces spurious logarithms that have no 
physical significance. Thus, although the MS mass is well-defined, it is not a particularly 
useful quantity to describe B decays. Therefore, several "threshold mass" definitions have 
been introduced that are more appropriate for low energy processes. 

The IS" mass is related to the pole mass through the perturbative relation [28, 29] 



7776 8 



,2 «3(/^) (, , 11 



(12) 



where the right hand side is the mass of the T(15') 66 bound state as computed in per- 
turbation theory, and I — ln[/x/(Q;5(/x) Cpmb)]. For the 15" mass there is a subtlety in the 
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perturbative expansion due to a mismatch between the order in e and the order in Ug, so 
that terms of order a^^^ in Eq. (12) are of order [28]. 

The potential-subtracted (PS) mass [30] is defined with respect to a factorization scale 
yUj. It is related to the pole mass through the perturbative relation 



m 



PS 



= 1 



rrib 



TT 



rrib 



le + €' 



27r 



+ y)/3o + ... 



(13) 



where now i = \n{fi/fif). In this paper we will choose fif = 2 GeV. 

Another popular definition is the kinetic, or "running", mass mb{fi) introduced in [25, 31]. 
The kinetic mass has properties similar to the PS mass, since it is defined with a cutoff that 
explicitly separates long- and short-distance physics. It should give comparable results, so 
we will not consider it here. We note, however, that in this scheme matrix elements such as 
Ai are also naturally defined with respect to a momentum cutoff. This has the advantage of 
absorbing some "universal" radiative corrections into the definitions of the matrix elements 
instead of the coefficients in the OPE, and is expected to improve the behavior of the 
perturbative series relating Ai to physical quantities. However, as usual, the perturbative 
relation between physical quantities is unchanged, and adopting this definition leaves our 
fits to \Vcb\ and nib unchanged. 

The results for the various shape variables are functions of the b quark mass. To simplify 
the expressions, in analogy with A defined in Eq. (5), we define new hadronic parameters 
by the following relations 
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A 



A^s = 4.2 GeV - mb(mb) . 



— m. 



15 



(14) 



We will refer to A, A^'^, A^^ and A^^^ generically as A. Note that the introduction of A is 
purely for computational convenience. The form Eq. (14) is chosen so that the value of A 
is numerically of order Aqcd- We can therefore expand the radiative corrections in powers 
of A and keep only the leading term and the first derivative. This is convenient because 
it avoids having to compute the radiative corrections, which involve a lengthy numerical 
integration, for each trial value of the quark mass in the fit. Note also that in the 15*, PS 
and MS schemes the dependence on — m;, is purely kinematic and is treated exactly, 
although it is formally of order Aqcd- 

Thus, the decay rates will be expressed in terms of 9 parameters: the A's in each mass 
scheme which we treat as order Aqcdj two parameters of order Aq^j-,, Ai and A2, and six 
parameters of order Aqqj-,, pi, p2, and 7i_4. Of these, only 6 are independent unknowns, as 
A2 is determined by Eq. (9), P2 — + %) is determined by Eq. (10), and T4 can be set to 
zero as explained preceding Eq. (7). 



IV. EXPANSIONS AND THEIR CONVERGENCE 

The computations in this paper include contributions of order I/ttlq and I/itlq, as well 
as radiative contribution of order e, and e^i^^, the so-called BLM contribution at order 
which is proportional to (3o- The dominant theoretical errors arise from the higher order 
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terms which we have neglected. In the perturbative series, we have neglected the non-BLM 
part of the two-loop correction. We have also neglected the unknown order fts/m^ and 
corrections in the OPE. The decay distributions depend on the charm quark mass, which 
is determined from ms — mo using Eq. (6). This formula introduces Aqc^/m^ corrections. 
Since rric only enters inclusive decay rates in the form m^/ml, the largest 1/m^ corrections 
are of order Aqq^/ {mlml). Finally, the O(eA) corrections for and S2 have only been 
calculated without a cut on the lepton energy [18]. 

For the B — > Xg^z/ decay rate and the shape variables defined in Eqs. (1), (2), and (4) we 
give results in the Appendix in the four different mass schemes discussed, for the coefficients 
X^^'^'^^Eq) in the expansion 



X{Eo) = X^'\Eo) + X^^\Eo) a + X^'\Eo) A^ + X^'\Eo) A^ 

+ X(5)(Eo) Ai + X^'\Eo) X2 + X^'\Eo) AiA + X^'\Eo) A2A 

+ X^'\E,) p, + X(i°)(i?o) P2 + X^''\Eo) T, + X^''\Eo) % ( 



where X(E^ is any of i?o(0,-Eo), Ri{Eo), S'j(-E'o)) or Tj(i?o) ^-nd i = 1,2. Note that to 
obtain Rq{Eq, Ei) one needs to reexpand i?o(0, i?i)/-Ro(0, Eq), but using -Ro(0, Eq) allows us 
to tabulate the results as a function of only one variable. The expressions for -Ro(0, Eq) are 
also convenient for deriving the predictions for other observables, such as those in Eq. (3). 

Unfortunately there is no simple way to relate the results in different mass schemes, 
because a particular value of the physical Eq cut corresponds to different limits of integrations 
in the dimensionless variables (such as 2EQ/mb) in different mass schemes. We list the 
coefficients of the expansions of the shape variables in the various mass schemes in the 
Appendix. 

Before using these expressions, one has to assess the convergence of both the perturbative 
expansions and of the power suppressed corrections. As each shape variable arises from a 
ratio of two series, the result can be worse or better behaved than the individual series in the 
numerator and denominator. We have checked that this is the reason for the apparent poor 
behavior of, for example, i?i(1.5 GeV) in the IS* scheme, where one sees that order ag term 
i?S^^^(1.5GcV) = 0.001, whereas the order ELM term i?S^^^(1.5 GeV) = 0.003 is larger. 
Since separately the numerator and denominator show good convergence, one should not 
conclude that i?i(1.5 GeV) is not a useful observable to constrain the HQET parameters. In 
general, one cannot conclude whether a series is poorly behaved or not by comparing the 
term with the term because of possible cancellations. Instead, one should compare with 
the expected size of terms based on a naive dimensional estimate. 

In Refs. [5, 18] the second hadronic invariant mass moment defined as ((m^ — m|,)^) 
was studied, and it was observed that the size of the Aq^q/ttt,^ correction was comparable 
to both the AQQj-,/m^ and ctsAQco/'^^b terms. The authors therefore concluded that the 
convergence of the OPE was suspect for this moment, and argued that useful constraints on 
A and Ai could not be obtained. A very similar situation holds for the variance ^2. However, 
one can obtain more insight into the convergence of this moment by examining the behavior 
of the relevant terms in the OPE for (m^) and {m%) separately. In the pole scheme (for 
simplicity), the expressions are 





+ 0.24 + 0.26^ + 1.02^ + 2.2-^ + 0.21-^ + 0.41 



as A 



E(>0 
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rrih A A Ai p^ a, a, A 



+ 0.07 + 0.14^ + 0.15^ - 0.23^ + 0.08-^ + 0.27 



Ee>o m% TUB ml ml ml ' 47r ' Anms 

(16) 



The OPE for both observables is well behaved, with the canonical size of the pi term a 
factor of 5-10 smaller than the Ai term. The corresponding constraints in the A — Ai plane 
have slopes which differ by roughly a factor of two, and so constrain one linear combination 
of A and Ai much better than the orthogonal combination. 

If instead of the second moment we consider the variance, we may combine the two series 
to find 

1 , /I , o V oi A^ Ai p^ a, a, A , , 
^(mi-(mi)2 =0.01^-0.14^-0.86-^ + 0.02-^ + 0.06-^ . (17 

The variance gives constraints in the A — Ai plane which are almost orthogonal to those of 
the first moment, but since it is simply a linear combination of the first and second moments, 
it cannot constrain the parameters any better. However, it is also no worse: none of the 
coefficients are larger than would be expected by dimensional analysis. The apparent poor 
convergence of the variance is due to a cancellation in the A (and to a lesser extent the Ai) 
terms between the two series. Therefore, there is no reason to expect the 0(l/m^) terms to 
be anomalously large. Constraints arising from 5*2 (or from {{m\ ~ ^i))^)) therefore need 
not be dismissed, although they are very sensitive to pi and so are of limited utility unless 
a sufficiently large number of observables is measured that pi is also constrained. 



V. EXPERIMENTAL DATA 

The experimental data for the lepton spectrum from the CLEO collaboration are the 
three lepton moments [9, 10] 

i?o(l-5GeV, 1.7GeV) = 0.6187 + 0.0021, 
i?i(1.5 GeV) = (1.7810 ± 0.0011) GeV, 

i?2(1.5 GeV) = (3.1968 ± 0.0026) GeV^ (18) 

For Rq and we used the averaged electron and muon values, with the full correlation 
matrix as given in Ref. [9]. For i?2, we have used the weighted average of the electron and 
muon data [10]. The DELPHI collaboration measures the lepton energy and variance [12], 

i?i(0) = (1.383 + 0.015) GeV, 

i?2(0) - Ri{Qf = (0.192 + 0.009) GeV^ (19) 

For the hadronic invariant mass spectrum we have CLEO measurements of the mean 
invariant mass and variance with a lepton energy cut of 1.5 GeV [8] 

,5i(1.5 GeV) = (0.251 + 0.066) GeV^ 

,^2(1.5 GeV) = (0.576 + 0.170) GeV^ (20) 

and DELPHI measurements of the mean invariant mass and variance with no lepton energy 
cut [13] 

5^(0) = (0.553 + 0.088) GeV^ 

^2(0) = (1.26 + 0.23) GeV^ (21) 
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Both collaborations also measure the second moment {{ni^ but we do not use this 

result since it is not independent of 5*1 and 5*2. 

The BABAR collaboration measures the first moment of the hadron spectrum for various 
values of the lepton energy cut [11]. The data points are highly correlated, and the variation 
of the first moment with the energy cut appears to be in poor agreement with the OPE 
predictions. We will do our fits without the BABAR data, as well as including the BABAR 
data for the two extreme values of their lepton energy cut, E = 0.9 and E = 1.5 GeV [11], 
to avoid overemphasizing many points with correlated errors in the fit, 

5i(1.5GeV) = (0.354 ± 0.080) GeV^ 

5i(0.9GeV) = (0.694 ± 0.114) GeVl (22) 

Note that we took into account that CLEG [9] and BABAR [11] used ttid = 1.975 GeV to 
obtain the quoted values of 5*1, whereas DELPHI [13] used mo = 1.97375 GeV. 
For the photon spectrum we use the GLEO results [7] 

ri(2GeV) = (2.346 ± 0.034) GeV, 

T2(2 GeV) = (0.0226 ± 0.0069) GeVl (23) 

The final piece of data is the semileptonic decay width, for which we use the average of 
S± and 5° data [32], 

T{B Xiu) = (42.7 ± 1.4) x 10"^^ MeV. (24) 

We do not average this with the Bg and 6-baryon semileptonic widths, as the power sup- 
pressed corrections can differ in these decays. 

Eqs. (18)-(24) provide a total of 14 measurements that enter our fit. 

YI. THE FIT 

In this section we perform a simultaneous fit to the various experimentally measured 
moments and the semileptonic rate. It is important to note that we do not include any 
correlations between experimental measurements beyond those presented in [9, 10], and so 
the experimental uncertainties are not completely taken into account. Nevertheless, the 
fit demonstrates the importance of including the full correlation of the 0(l/mf) terms in 
the different observables, and also indicates the relative importance of the theoretical and 
experimental uncertainties. 

We use the fitting routine Minuit to fit simultaneously for the shape variables and the 
total semileptonic branching fraction, by minimizing x^, and present results for the fit in 
the IS scheme (the other schemes give comparable results). 

In addition to the experimental imcertainties, there are also uncertainties in the theory 
because the formulae used in the fit are not exact. From naive dimensional analysis we 
find the fractional theory errors 0.0003 from (as/^ir)'^ terms, 0.0002 from {as/4:7r) Aqcc/^fe 
terms, and 0.001 from Aqqj-, /(m^m^) terms. In some cases, naive dimensional analysis 
underestimates the uncertainties, and an alternative is estimating the uncertainties by the 
size of the last term computed in the perturbation series. We combine these estimates by 
adding in quadrature half of the e^Lj^ term and a O.OOlm^ theoretical error for quantities 
with mass dimension n. In computing x^, we add this theoretical error in quadrature to the 
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TABLE I: Fit results for \Vcb\, mb, Ai and Ai + (71 + 372)/mfe in the l^* scheme. The iV^cfel 
value includes electromagnetic radiative corrections; see Eq. (26). The upper/lower blocks are fits 
excluding/including the BABAR data, and have 5 and 7 degrees of freedom, respectively. 



[GeV] 


x' 


\Vcb\ X 10^ 


mf [GeV] 


Ai [GeV^] 


Ai + ^ [GeV^] 


0.5 
1.0 


5.0 
3.5 


40.8 ± 0.9 
41.1 ±0.9 


4.74 ±0.10 
4.74 ±0.11 


-0.22 ±0.38 
-0.40 ±0.26 


-0.31 ±0.17 
-0.31 ±0.22 


0.5 
1.0 


12.9 
8.5 


40.8 ± 0.7 

40.9 ± 0.8 


4.74 ±0.10 
4.76 ±0.11 


-0.14 ±0.13 
-0.22 ±0.25 


-0.29 ±0.10 
-0.17 ±0.21 



TABLE II: Fit results for the 1/m^ coefficients in the IS" scheme. The upper /lower blocks are fits 
excluding/including the BABAR data. The constraint in Eq. (10) is used to determine p2- 



[GeV] 


Pi 


P2 


T1+T3 


Ti + 3T2 


0.5 


0.15 ±0.12 


-0.01 ±0.11 


-0.15 ±0.84 


-0.45 ± 1.11 


1.0 


0.16 ±0.18 


-0.05 ± 0.16 


0.11 ± 0.10 


0.-15 ± 0.49 


0.5 


0.17 ±0.09 


-0.04 ± 0.09 


-0.34 ±0.16 


-0.66 ± 0.32 


1.0 


0.08 ±0.18 


-0.12 ±0.15 


0.11 ±0.33 


0.23 ± 0.47 



experimental errors. This procedure avoids giving a large weight in the fit to a very accurate 
measurement that cannot be computed reliably. Because the perturbative results in the l^S" 
scheme are not expected to be artificially badly behaved (as they are in, for example, the 
pole scheme) this estimate of the perturbative uncertainty should be reasonable. We will 
examine the convergence of perturbation theory later in this section. 

The unknown matrix elements of the l/ni'^ operators are the largest source of uncertainty 
in the fit. One expects these matrix elements to be of order Aq^q. To allow for this 
theoretical input, we include an additional contribution to from the matrix elements of 
each 1/ml operators, pi_2 and Ti_4, that we denote generically by {O), 

fO, \{0)\<ml, 

where {m^^, M^) are both thought of as quantities of order Aqcd- This way we do not 
prejudice {O) to have any particular value in the range 1(C) I < ml- In the fit we take 
= 500 MeV, and vary between 500 MeV and 1 GeV to test that our results for \Vcb\ 
and rrib are insensitive to this input (our final results are obtained with = 500 MeV). The 
data are sufficient to constrain the 1/ml operators in the sense that they can be consistently 
fit with reasonable values, but they are not determined with any useful precision. Finally, 
since only three linear combinations of 7^_4 appear in the formulae, we fit setting T4 — 0, 
so that the fit values for 7^_3 with this choice for T4 are the values of 7^ — 37^, T2 + T4, and 
T3 ± 3%. 

The fit results are summarized in Tables I and IL In Table I we show the results of the fit 
for \Vcb\, and Ai , as well as the "effective" combination Ai + (71 + 37^) / rub which enters 
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FIG. 1: Comparison of the BABAR measurement of the hadron invariant mass spectrum [11] 
vs. the lepton energy cut (black squares), and our prediction from the fit not including BABAR 
hadronic mass data (red triangles). 



in the OPE, and which, due to correlated errors, is better constrained than Ai. From these 
results we can also obtain an expression for |V"cb| as a function of the semileptonic branching 
ratio and the B meson lifetime. We find 



\Vcb\ = (40.8 ±0.7) X lO-^r^QED 



B{B ^XJu) 1.6ps 



.1/2 



0.105 



(26) 



The quoted error contains all uncertainties from m^, Ai, the 1/m^ matrix elements, as well as 
perturbative uncertainties. The parameter i/qed ~ 1.007 is the electromagnetic correction 
to the inclusive decay rate, which has been included in the values for |Vcfe| presented in 
Table I. Including the BABAR data increases the by about a factor of two. Doubling 
the allowed range of the 1/mf parameters increases the uncertainties only minimally and 
reduces somewhat. 

The reason we carried out separate fits excluding and including the BABAR data on 
Si{Eq) is because of its inconsistency at low Eq with the fit done without it. To see this, 
note that on very general grounds 5*1(^^0) is a monotonically decreasing function of £"0- 
The theoretical prediction corresponding to the fit in the first line of Table I is 5*1(0) = 
(0.42±0.03) GeV^ which is significantly below the lowest BABAR data point, ^i(0.9 GeV) = 
(0.694 ± 0.114) GeV. Assuming that the branching ratio to nonresonant channels between 
D* and D** is negligible, this prediction for S'i(O) implies an upper bound on the fraction 
of excited (i.e., non-D^*^) states in S — ^ Xciu decay [16], which is below 25%, and is in 
contradiction with the measured B D'^*Hi' branching fractions. To resolve this, either 
the assumption that low-mass nonresonant channels are negligible could be wrong, or some 
measurements or the theory have to be several standard deviations off. The X^. spectrum 
effectively has this (assumed) feature in the CLEO and BABAR analyses but not in DELPHI. 
It is thus crucial to precisely and model independently measure the mx^ distribution in 
semileptonic B —>■ Xc^v decay. A comparison of the BABAR hadronic moment data with 
our fit is given in Fig. 1. 
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TABLE III: Fit predictions for fractional moments of the electron spectrum. The upper/lower 
blocks are fits excluding/including the BABAR data. 



[GeV] 




R3b 




R^b 


D3 


D4 


0.5 


0.302 ± 0.003 


2.261 ± 0.013 


2.127 ±0.013 


0.684 ± 0.002 


0.520 ± 0.002 


0.604 ± 0.002 


1.0 


0.302 ± 0.002 


2.261 ±0.011 


2.128 ±0.011 


0.684 ± 0.002 


0.519 ±0.002 


0.604 ± 0.001 


0.5 


0.302 ± 0.002 


2.261 ±0.012 


2.127 ±0.012 


0.684 ± 0.002 


0.520 ± 0.002 


0.604 ± 0.001 


1.0 


0.302 ± 0.002 


2.262 ± 0.012 


2.129 ±0.012 


0.684 ± 0.002 


0.519 ±0.001 


0.604 ± 0.001 



To get more insight into the obtained uncertainties, we have performed several additional 
fits in which we turn off individual contributions to the errors. Here we present the results 
for the fits with — 0.5 and not including the BABAR data. Similar results are true when 
the BABAR values are included. Neglecting all 1 / mf terms, as well as the naive estimate of 
the theoretical uncertainties gives a fit with = 81 for 9 degrees of freedom. Including only 
the 1/ml terms gives = 21 for 5 degrees of freedom. This is a vastly better fit, reducing 
by about 60 by adding only 4 new parameters. Nevertheless, the fact that per degree of 
freedom is about 5 shows that there is a statistically significant discrepancy between theory 
and experiment if other theoretical uncertainties are not included. Only after including this 
estimate do we get x^/dof ~ 1. We also estimated the size of the theoretical uncertainties 
by setting all experimental errors to zero. This reduces all uncertainties by roughly a factor 
of three. Thus, the fit is dominated by experimental uncertainties. 

The fit gives a value of the b quark mass which is consistent with other extractions, 
and with an uncertainty at the 100 MeV level. For comparison, T sum rules extractions in 
Refs. [33, 34] give mf = 4.69 ± 0.03 GeV and m^^ = 4.78 ± 0.11 GeV, respectively by a 
fit to the BB system near threshold. The error on Ai is larger than previous extractions 
from Ti and 5*1 [8], because we arc inchiding more conservative estimates of the theoretical 
uncertainties. Despite this, the uncertainty on \ Vch\ is smaller than from previous extractions. 
Note that we have only used the value of the semileptonic branching ratio of B mesons. It 
is inconsistent to combine the average semileptonic branching ratio of b quarks (including 
Bg and A;, states) with the moment analyses, since hadronic matrix elements have different 
values in the B/B* system, and in the Bg/B* or A;,. 

The fit results for the 1/ml parameters are shown in Table II. Clearly, one is not able 
to determine the values of the 1/ml parameters from the present fit. All that can be said 
is that the preferred values are consistent with dimensional estimates. There is also some 
indication that p2 is small, as is expected in some models [16]. 

One can also use the fits to predict other observablcs that can be measured. For example, 
we predict the values for the fractional moments R^a, Rsb, R^a-i R'lb, -D3 and D4 given by 
Bauer and Trott [35]. The predicted values are given in Table III. The results are robust, 
and do not depend on the width chosen for the operators, or whether or not we include 
the BABAR data. 

Finally, it is useful to study the convergence of perturbation theory by carrying out the 
fit at different orders in the perturbation expansion. In Figure 2 we show the la error ellipse 
in the ml^ vs. \Vcb\ plane, for the four different mass schemes. For each scheme we show 
three contours, obtained at tree level (dotted red curves), at order e (dashed blue curves). 
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FIG. 2: The la error ellipse in the ml vs. \Vcb\ plane, using different mass schemes for the fit. For 
each scheme we show the contours obtained at tree level (dotted red curves), at order e (dashed 
blue curves), and at order e^Ljyj (solid black curves). 

and including order e^L^ corrections as well (solid black curves). For each of these curves, 
the conversion of the fitted mass to the IS" mass has been done at the consistent order 
in perturbation theory. One can see that the convergence of the perturbative expansion 
is slightly better for the 15* and the PS schemes compared with the pole scheme. This is 
because there is an incomplete cancellation of formally higher order terms, such as a^A^, 
which are large in the pole scheme. The larger uncertainties in the MS scheme are due 
to large contributions at BLM order, which are included in the uncertainty estimate, as 
explained at the beginning of this section. 



VII. SUMMARY AND CONCLUSIONS 

Experimental studies of the shape variables discussed in this paper are crucial in deter- 
mining from experimental data the accuracy of the theoretical predictions for inclusive B 
decays rates, which rest on the assumption of local duality. Detailed knowledge of how well 
the OPE works in different regions of phase space (and a precise value of mi,) will also be 
important for the determination of \Vub\ from inclusive B decays. A serious discrepancy 
between theory and data would imply, for example for |V"cf,|, that only its determination 
from exclusive decays has a chance of attaining a reliable error below the ~ 5% level. 

The analysis in this paper shows that at the present level of accuracy, the data from 
the lepton and photon spectra are consistent with the theory, with no evidence for any 
breakdown of quark-hadron duality in shape variables. Two related problems at present 
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FIG. 3: The 1- and 2-a regions in the vs. \Vcb\ plane using the IS* mass scheme. Superimposed 
are the values and errors of the determination of \Vcb\ from exclusive decays [32] and that of ml^ 
from sum rules in Ref. [33] (red square) and Ref. [34] (magenta triangle). 



are the BABAR measurement of the average hadronic invariant mass as a function of the 
lepton energy cut and the total branching fraction to D and D* states, both of which appear 
problematic to reconcile with the other measurements combined with the OPE. However, 
both problems depend on assumptions about the invariant mass distribution of the decay 
products, which needs to be better understood. Excluding the BABAR data and the problem 
of the B D^*^^ branching ratios, the fit provides a good description of the experimental 
results, with = 5.0 for 12 data points and 7 fit parameters in the 15* scheme. 

The main results (in the IS* scheme) are summarized in Fig. 3 where we compare our 
determination of \Vcb\ and ml^ with those from exclusive B decays and upsilon sum rules. 
We obtain the following values: 



IK 



m 



cb\ 
IS 



(40.8 ±0.9) X 10"^ 
(4.74 ±0.10) GeV. 



(27) 



This corresponds to the MS mass mb{mb) = 4.22 ± 0.09 GeV. We have also presented the 
value of \Vcb\ as a function of the semileptonic branching ratio and the B meson lifetime 



\Vcb\ = (41.1 ±0.7) X 10' 



B{B XJi?) 1.6ps 



.1/2 



0.105 



(28) 



We have constrained the 1/m^ matrix elements and predicted the values for fractional mo- 
ments of the electron spectrum to better than 1% accuracy. 

Setting experimental errors to zero gives errors in \Vcb\ and ml^ of 0.35 x 10~^ and 
35 MeV, respectively. These numbers indicate the theoretical limitations, although their 
precise values depend on details of how the theoretical uncertainties are estimated. If the 
agreement between the experimental results improve in the future, then a full two loop 
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calculation of the total semileptonic rate and of 5 — > Xc£u decay spectra would help to 
further reduce the theoretical uncertainty in \ Vcb\ and mi,. 
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APPENDIX: COEFFICIENT FUNCTIONS IN VARIOUS MASS SCHEMES 

In this Appendix we give numerical results for the the B —>■ XcIp decay rate and the 
shape variables defined in Eqs. (1), (2), and (4), in the four mass schemes discussed. For 
all quantities the coefficients of the expansions are defined as in Eq. (15), and all numerical 
values are in units of GeV to the appropriate power. We use as{mb) = 0.22 and the spin- 
and isospin-averaged meson masses, = 5.314 GeV and mu = 1.973 GeV. 



1. The \S mass scheme 



The B — > Xc(-v decay width in the 15" scheme is given by 



T{B^ XJP) = 



Gl\V^fmr\' 



0.534 - 0.232 A - 0.023 A^ + 0. A^ 



1927r3 V 2 y 
-0.11 Ai - 0.15 A2 - 0.02 AiA + 0.05 A2A 

-0.02 pi + 0.03 p2 - 0.05 Ti + 0.01 Ta 

-0.07 % - 0.03 % - 0.051 e - 0.016 cIlm + 0.016 eA 



(A.l) 



We tabulate the shape variables defined in Eq. (1) in Tables IV, V, and VI, and those 
defined in Eq. (2) in Tables VII and VIII in the IS* mass scheme. For and S2 we do not 
show the £^o-dependence of the order eA terms, as they are not known. For all quantities 
the coefficients of the expansions are defined as in Eq. (15). 

For the B X,^ shape variables defined in Eq. (4), only tI^^\ tI;^^\ and T^^^^ are 
functions of Eq, once 771^/2 — Aqcd- For the other T"s in the 15" scheme we find 



n 



(1) 



mr 



tP = -0.01, 



n 
-0.01, 



(3) 



(8) 



1 

2' 

-0.03, Tl 



.(9) 



(4) 



0, Tf ^ = -0.05, TP = -0.16, 



(6) 



-0.02, 
-0.03, 



(10) 



0.18, 



(A.2) 
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TABLE IV: Coefficients for Rq{0, Eq) in the 15" scheme as a function of Eq. 



Eq 




Kg 


Kg 


Kg 


Kg 


d(6) 

Kg 


Kg 


Kg 


Kg 


d(10) 

Kg 


Kg 


Kg 


d(13) 

Kg 


Kg 


d(1B) 

Kg 


d(16) 

Kg 


Kg 


0.5 


0.972 


-0.003 


-0.002 


0. 


0. 


-0.01 


0. 


0. 


0. 


0. 


0. 


0. 





0. 


0. 


0. 


0. 


0.7 


0.927 


-0.008 


-0.005 


0. 


-0.01 


-0.03 


-0.01 


-0.01 


0. 


0. 


-0.01 


0. 


-0.01 


-0.01 


0.001 


0.001 


0. 


0.9 


0.853 


-0.016 


-0.01 


-0.01 


-0.02 


-0.06 


-0.02 


-0.03 


0. 


0.01 


-0.01 


0. 


-0.02 


-0.01 


0.002 


0.001 


0. 


1.1 


0.749 


-0.028 


-0.015 


-0.01 


-0.04 


-0.1 


-0.03 


-0.05 


-0.01 


0.01 


-0.02 


0. 


-0.03 


-0.02 


0.002 


0.001 


0. 


1.3 


0.615 


-0.043 


-0.022 


-0.01 


-0.06 


-0.15 


-0.05 


-0.08 


-0.01 


0.02 


-0.03 


0. 


-0.04 


-0.03 


0.003 


0.002 


0. 


1.5 


0.455 


-0.062 


-0.029 


-0.01 


-0.08 


-0.2 


-0.07 


-0.11 


-0.01 


0.03 


-0.04 


0. 


-0.05 


-0.04 


0.003 


0.002 


0. 


1.7 


0.279 


-0.084 


-0.037 


-0.02 


-0.1 


-0.25 


-0.08 


-0.15 


-0.01 


0.03 


-0.04 


-0.01 


-0.06 


-0.05 


0.002 


0.003 


-0.001 



TABLE V: Coefficients for Ri{Eq) in the IS scheme as a function of Eq. 



Ei, 
















Rf^ 


Rf^ 






r['^^ 


r[''^ 


r[''^ 


r['"^ 




r[''^ 





1.392 


-0.077 


-0.026 


-0.01 


-0.11 


-0.22 


-0.07 


-0.08 


-0.04 


0.01 


-0.04 


-0.02 


-0.05 


-0.05 


0.003 


0.003 


0. 


0.5 


1.422 


-0.076 


-0.025 


-0.01 


-0.11 


-0.22 


-0.06 


-0.08 


-0.04 


0.01 


-0.04 


-0.02 


-0.05 


-0.05 


0.003 


0.003 


0. 


0.7 


1.461 


-0.075 


-0.023 


-0.01 


-0.11 


-0.21 


-0.06 


-0.08 


-0.04 


0.01 


-0.04 


-0.02 


-0.05 


-0.04 


0.002 


0.003 


0. 


0.9 


1.517 


-0.074 


-0.022 


-0.01 


-0.11 


-0.2 


-0.06 


-0.08 


-0.04 


0.01 


-0.04 


-0.02 


-0.05 


-0.04 


0.002 


0.003 


0. 


1.1 


1.588 


-0.074 


-0.021 


-0.01 


-0.11 


-0.19 


-0.06 


-0.08 


-0.04 


0. 


-0.04 


-0.03 


-0.04 


-0.04 


0.001 


0.003 


0. 


1.3 


1.672 


-0.075 


-0.02 


-0.01 


-0.11 


-0.19 


-0.06 


-0.07 


-0.05 


0. 


-0.04 


-0.03 


-0.04 


-0.04 


0.001 


0.003 


0. 


1.5 


1.767 


-0.077 


-0.02 


-0.01 


-0.12 


-0.17 


-0.07 


-0.07 


-0.06 


-0.02 


-0.04 


-0.04 


-0.04 


-0.04 


0.001 


0.003 


0. 


1.7 


1.872 


-0.08 


-0.021 


-0.01 


-0.14 


-0.16 


-0.1 


-0.06 


-0.1 


-0.04 


-0.04 


-0.06 


-0.03 


-0.03 


0.001 


0.003 


0. 



TABLE VL Coefficients for R2{Eq) in the 15 scheme as a function of Eq. 



Eq 


<^ 


o(2) 
K^ 


o(3) 


o(4) 

rtj 


p(5) 


p(6) 


p(7) 


p(8) 
K^ 


p(9) 

K2 


p(10) 

K2 


r(11) 


r(12) 


p(13) 
K2 


K2 


p(16) 
K2 


p(16) 
K2 


p(17) 
K2 





2.118 


-0.247 


-0.07 


-0.02 


-0.36 


-0.68 


-0.19 


-0.21 


-0.15 


0.02 


-0.14 


-0.08 


-0.16 


-0.14 


0.008 


0.01 


-0.001 


0.5 


2.175 


-0.247 


-0.069 


-0.02 


-0.36 


-0.68 


-0.19 


-0.22 


-0.15 


0.02 


-0.13 


-0.08 


-0.16 


-0.14 


0.007 


0.009 


-0.001 


0.7 


2.263 


-0.248 


-0.067 


-0.02 


-0.36 


-0.68 


-0.19 


-0.22 


-0.16 


0.01 


-0.13 


-0.09 


-0.15 


-0.14 


0.007 


0.009 


-0.001 


0.9 


2.401 


-0.252 


-0.065 


-0.02 


-0.37 


-0.67 


-0.19 


-0.22 


-0.17 


0.01 


-0.13 


-0.1 


-0.15 


-0.14 


0.005 


0.009 


-0.001 


1.1 


2.593 


-0.259 


-0.064 


-0.02 


-0.38 


-0.67 


-0.19 


-0.23 


-0.18 


-0.01 


-0.13 


-0.11 


-0.14 


-0.14 


0.004 


0.009 


-0.001 


1.3 


2.842 


-0.271 


-0.063 


-0.02 


-0.41 


-0.66 


-0.21 


-0.23 


-0.21 


-0.03 


-0.14 


-0.13 


-0.14 


-0.14 


0.003 


0.009 


-0.001 


1.5 


3.15 


-0.288 


-0.066 


-0.02 


-0.46 


-0.64 


-0.24 


-0.23 


-0.28 


-0.07 


-0.14 


-0.18 


-0.13 


-0.14 


0.003 


0.011 


-0.001 


1.7 


3.518 


-0.311 


-0.072 


-0.02 


-0.58 


-0.62 


-0.35 


-0.21 


-0.43 


-0.16 


-0.16 


-0.26 


-0.12 


-0.13 


0.004 


0.013 


0. 



and 

^(1) ^ rp{2) ^ rp{3) ^ ^(4) ^ ^(6) ^ y(7) ^ ^(8) ^ ^(13) ^ ^(14) ^ q 

= ^ = -0.04, = -T^^^^ = 0.05, t^^^^ = -0.02. 

The remaining, £^o-dependent coefficients of the perturbative corrections are fisted in Ta- 
ble IX. 



2. The PS mass scheme 

The expressions for the B XJm decay rate and the shape variables in the PS scheme 
are almost identical to Eq. (A.l), Tables IV-VIII, and Eqs. (A. 2) and (A. 3), because we 
choose to expand m^^ about m-^jl as well. The difference in the B — >■ Xf.iv rate compared 
with Eq. (A.l) is that the perturbation series is replaced by —0.020 e — 0.003 cIlm 

+ 0.025 eA, 

and of course, the meaning of A changes from A^'^ to A^^. 
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TABLE VII: Coefficients for Si{Eo) in the l^* scheme as a function of Eq. 





s[' 


s[' 




s[' 






sY' 


S^' 


SY> 


si"' 


S\ > 


Si > 




s[ > 


s[^' 




s["' 





0.832 


1.633 


0.416 


0.13 


1.49 


-0.36 


0.75 


0. 


0.46 


-0.24 


0.53 


0.25 


0.5 


0.14 


0.044 


-0.025 


0.025 


0.5 


0.82 


1.609 


0.409 


0.12 


1.5 


-0.32 


0.75 


0.02 


0.48 


-0.24 


0.54 


0.26 


0.5 


0.14 


0.039 


-0.028 




0.7 


0.805 


1.578 


0.398 


0.12 


1.52 


-0.26 


0.77 


0.05 


0.5 


-0.23 


0.54 


0.27 


0.5 


0.16 


0.032 


-0.031 




0.9 


0.784 


1.533 


0.38 


0.11 


1.56 


-0.16 


0.79 


0.12 


0.55 


-0.22 


0.55 


0.3 


0.51 


0.18 


0.023 


-0.035 




1.1 


0.759 


1.479 


0.354 


0.1 


1.63 


-0.02 


0.83 


0.22 


0.63 


-0.2 


0.57 


0.34 


0.52 


0.2 


0.011 


-0.04 




1.3 


0.734 


1.42 


0.319 


0.09 


1.74 


0.18 


0.91 


0.38 


0.77 


-0.16 


0.59 


0.41 


0.54 


0.24 


-0.002 


-0.046 




1.5 


0.716 


1.371 


0.277 


0.06 


1.97 


0.45 


1.07 


0.65 


1.03 


-0.06 


0.64 


0.55 


0.56 


0.3 


-0.018 


-0.054 




1.7 


0.72 


1.368 


0.254 


0.05 


2.49 


0.84 


1.59 


1.13 


1.64 


0.22 


0.76 


0.86 


0.6 


0.38 


-0.035 


-0.066 





TABLE VIII: Coefficients for S2{Eq) in the IS* scheme as a function of Eq. 







o(2) 
•^2 


A3) 
■^2 






•^2 


o(7) 
•^2 




o(9) 
•^2 


AW) 

■^2 


■^2 








•^2 


o(16) 
■^2 







0.125 


0.472 


0.531 


0.16 


-4.43 


-0.68 


-1.04 


-1.6 


-5.46 


1.07 


-0.94 


-2.8 


-0.05 


-0.13 


0.381 


-0.428 


0.171 


0.5 


0.123 


0.467 


0.524 


0.16 


-4.34 


-0.66 


-0.99 


-1.55 


-5.53 


0.96 


-0.93 


-2.74 


-0.05 


-0.12 


0.405 


-0.42 




0.7 


0.123 


0.465 


0.521 


0.16 


-4.23 


-0.64 


-0.91 


-1.5 


-5.64 


0.81 


-0.9 


-2.67 


-0.05 


-0.12 


0.448 


-0.408 




0.9 


0.124 0.468 


0.524 


0.16 


-4.08 


-0.62 


-0.78 


-1.43 


-5.85 


0.59 


-0.87 


-2.58 


-0.05 


-0.11 


0.526 


-0.391 




1.1 


0.126 


0.477 


0.533 


0.16 


-3.89 


-0.6 


-0.6 


-1.36 


-6.2 


0.28 


-0.83 


-2.46 


-0.05 


-0.11 


0.661 


-0.37 




1.3 


0.128 


0.486 


0.546 


0.17 


-3.69 


-0.57 


-0.35 


-1.28 


-6.79 


-0.11 


-0.79 


-2.33 


-0.05 


-0.1 


0.892 


-0.344 




1.5 


0.128 


0.487 


0.55 


0.18 


-3.5 


-0.53 


-0.04 


-1.19 


-7.88 


-0.61 


-0.75 


-2.21 


-0.05 


-0.1 


1.328 


-0.311 




1.7 


0.12 


0.454 0.509 


0.16 


-3.46 


-0.49 


0.16 


-1.08 


-10.34 


-1.34 


-0.74 


-2.18 


-0.05 


-0.09 


2.345 


-0.273 





TABLE IX: Perturbative coefficients for Ti{Eq) and T2{Eq) in the IS scheme as a function of Eq. 



Eo 


rp(15) 
^1 


7^(16) 
^ 1 


^(17) 


^2 


7.(16) 
^2 




1.7 


-0.043 


-0.017 


0.016 


0.016 


0.011 


-0.014 


1.8 


-0.038 


-0.014 


0.021 


0.012 


0.009 


-0.014 


1.9 


-0.032 


-0.011 


0.026 


0.01 


0.007 


-0.014 


2 


-0.025 


-0.006 


0.033 


0.007 


0.006 


-0.013 


2.1 


-0.017 


-0.001 


0.042 


0.004 


0.004 


-0.012 


2.2 


-0.007 


0.008 


0.056 


0.002 


0.002 


-0.01 



Next wc tabulate the coefficients of the perturbation series of the shape variables defined 
in Eqs. (1) and (2), that differ from the entries in Tables IV-VIII, in Table X in the PS 
mass scheme. For 5'i and S2 we do not show in the tables the order eA terms again as their 
£^o-dependence is not known. For all quantities the coefficients of the expansions are defined 
as in Eq. (15). 

For the B — > X57 shape variables defined in Eq. (4), the expressions for T2 are identical in 
the 15" and PS schemes, and so only Ti^''\ Ti^^\ and t/^^^ differ between these two schemes. 
The results for these coefficients in the PS scheme are shown in Table XI. 
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TABLE X: Perturbative coefficients for Ro{0,Eo), Ri{Eo), i?2(-Eo), Si{Eo), and S2iEo) in tlie PS 
scheme, that differ from the results in the IS scheme, as a function of Eq. 



Eo 










R\ ' 


R\ ' 








s[ ' 


Si ' 


^(17) 
SI ' 


sr' 


sr' 













0.013 


0.007 


0.008 


0.041 


0.024 


0.021 


-0.178 


-0.106 


-0.106 


0.317 


-0.452 


0.022 


0.5 


0.001 


0. 


0.001 


0.013 


0.007 


0.007 


0.041 


0.023 


0.02 


-0.18 


-0.108 




0.342 


-0.443 




0.7 


0.002 


0.001 


0.002 


0.012 


0.007 


0.007 


0.04 


0.023 


0.02 


-0.182 


-0.109 




0.385 


-0.432 




0.9 


0.004 


0.002 


0.003 


0.012 


0.007 


0.007 


0.04 


0.023 


0.019 


-0.186 


-0.111 




0.462 


-0.415 




1.1 


0.006 


0.003 


0.005 


0.011 


0.007 


0.006 


0.039 


0.024 


0.019 


-0.19 


-0.114 




0.596 


-0.393 




1.3 


0.009 


0.004 


0.007 


0.011 


0.007 


0.006 


0.04 


0.025 


0.019 


-0.195 


-0.116 




0.826 


-0.368 




1.5 


0.011 


0.006 


0.009 


0.011 


0.007 


0.006 


0.042 


0.027 


0.02 


-0.205 


-0.122 




1.262 


-0.335 




1.7 


0.013 


0.007 


0.01 


0.012 


0.008 


0.007 


0.046 


0.031 


0.023 


-0.221 


-0.135 




2.283 


-0.296 





TABLE XL Perturbative coefficients for Ti{Eo) in the PS scheme as a function of Eq. 











1.7 


0.025 


0.011 


0.022 


1.8 


0.03 


0.014 


0.026 


1.9 


0.036 


0.018 


0.032 


2 


0.043 


0.022 


0.038 


2.1 


0.051 


0.028 


0.047 


2.2 


0.061 


0.036 


0.062 



TABLE XIL Coefficients for ii!o(0, ^'o) in the MS scheme as a function of Eq. 



Eo 


o(l) 
Kg 


o(2) 

Kg 


o(3) 
Kg 


p(4) 


p(5) 


r(6) 

Kg 


p(7) 

Kg 


p(8) 

Kg 


p(9) 

Kg 


p(10) 

Kg 


r(11) 

Kg 


p(12) 


p(13) 

Kg 


p(14) 

Kg 


Kg 


p(16) 

Kg 


Kg 


0.5 


0.969 


-0.007 


-0.005 


0. 


-0.01 


-0.01 


-0.01 


-0.01 


0. 


0. 


-0.01 


0. 


-0.01 


0. 


0.003 


0.001 


0.003 


0.7 


0.92 


-0.017 


-0.013 


-0.01 


-0.02 


-0.04 


-0.02 


-0.02 


0. 


0.01 


-0.01 


0. 


-0.02 


-0.01 


0.007 


0.004 


0.008 


0.9 


0.841 


-0.033 


-0.025 


-0.02 


-0.04 


-0.07 


-0.04 


-0.05 


0. 


0.01 


-0.03 


0. 


-0.03 


-0.02 


0.013 


0.007 


0.015 


1.1 


0.729 


-0.054 


-0.04 


-0.03 


-0.06 


-0.13 


-0.07 


-0.09 


0. 


0.02 


-0.05 


0. 


-0.06 


-0.03 


0.021 


0.012 


0.023 


1.3 


0.584 


-0.08 


-0.056 


-0.04 


-0.1 


-0.2 


-0.11 


-0.14 


0. 


0.03 


-0.07 


0. 


-0.08 


-0.05 


0.031 


0.017 


0.032 


1.5 


0.411 


-0.11 


-0.071 


-0.05 


-0.13 


-0.29 


-0.15 


-0.22 


0. 


0.04 


-0.09 


0. 


-0.11 


-0.07 


0.041 


0.024 


0.039 


1.7 


0.221 


-0.145 


-0.086 


-0.05 


-0.16 


-0.36 


-0.18 


-0.3 


0.01 


0.04 


-0.11 


-0.01 


-0.13 


-0.09 


0.052 


0.032 


0.046 



3. The MS mass scheme 



The B — > Xciiy decay width in the MS scheme is given by 

|2 



r(B^ xjp) = 



G% \ Vcb\ 



(4.2 GeV)' 



1927r3 

-0.22 Ai - 0.22 A2 - 0.04 AiA + 0.1 A2A 
-0.01 pi + 0.05 p2 - 0.16 Ti + 0.01 r2 



0.733 - 0.464 A - 0.036 A^ + 0.01 A'^ 



-0.18 % - 0.05 T4 + 0.085 e + 0.065 e 



BLM 



0.022 eA 



(A.4) 



We tabulate the shape variables defined in Eq. (1) in Tables XII, XIII, and XIV, and 
those defined in Eq. (2) in Tables XV and XVI in the MS mass scheme. For 5*1 and 5*2 
we do not show the i^o-dependcncc of the order eA terms, as they are not known. For all 
quantities the coefficients of the expansions are defined as in Eq. (15). 

For the B — > Xg'y shape variables defined in Eq. (4), tI;^\ . . . Tj;^^^ are independent of Eq, 
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TABLE XIII: Coefficients for Ri{Eq) in the MS seheme as a function of Eq. 





















R'f^ 






Ri'''^ 










r["^ 





1.342 


-0.117 


-0.054 


-0.03 


-0.16 


-0.27 


-0.12 


-0.12 


-0.03 


0.01 


-0.09 


-0.03 


-0.1 


-0.07 


0.043 


0.026 


0.027 


0.5 


1.373 


-0.113 


-0.05 


-0.03 


-0.15 


-0.27 


-0.12 


-0.12 


-0.04 


0.01 


-0.09 


-0.03 


-0.1 


-0.06 


0.042 


0.025 


0.025 


0.7 


1.413 


-0.11 


-0.047 


-0.02 


-0.15 


-0.26 


-0.11 


-0.12 


-0.04 


0.01 


-0.09 


-0.03 


-0.1 


-0.06 


0.04 


0.024 


0.023 


0.9 


1.47 


-0.106 


-0.043 


-0.02 


-0.15 


-0.26 


-0.11 


-0.12 


-0.04 


0.01 


-0.08 


-0.03 


-0.09 


-0.06 


0.039 


0.024 


0.021 


1.1 


1.542 


-0.104 


-0.039 


-0.02 


-0.15 


-0.25 


-0.11 


-0.12 


-0.05 


0. 


-0.08 


-0.04 


-0.09 


-0.06 


0.037 


0.023 


0.019 


1.3 


1.626 


-0.103 


-0.036 


-0.02 


-0.15 


-0.24 


-0.11 


-0.12 


-0.06 


-0.01 


-0.08 


-0.05 


-0.08 


-0.06 


0.037 


0.023 


0.017 


1.5 


1.72 


-0.105 


-0.035 


-0.01 


-0.17 


-0.22 


-0.13 


-0.12 


-0.08 


-0.03 


-0.08 


-0.07 


-0.07 


-0.05 


0.037 


0.024 


0.016 


1.7 


1.823 


-0.109 


-0.036 


-0.01 


-0.22 


-0.2 


-0.22 


-0.1 


-0.16 


-0.08 


-0.08 


-0.11 


-0.06 


-0.05 


0.039 


0.025 


0.017 



TABLE XIV: Coefficients for R2{Eq) in the MS scheme as a function of Eq. 



Eo 


r,(l) 
R2 


R2 


r,(3) 
R2 


R2 


r(5) 


R2 


o(7) 
R2 


R2 


R2 


R2 


R2 


R2 


R2 




R2 


p(16) 
R2 







1.963 


-0.35 


-0.136 


-0.07 


-0.49 


-0.82 


-0.31 


-0.31 


-0.14 


0.02 


-0.27 


-0.11 


-0.3 


-0.2 


0.129 


0.077 


0.063 


0.5 


2.02 


-0.348 


-0.132 


-0.06 


-0.49 


-0.82 


-0.31 


-0.31 


-0.15 


0.02 


-0.27 


-0.11 


-0.3 


-0.2 


0.127 


0.077 


0.061 


0.7 


2.108 


-0.346 


-0.127 


-0.06 


-0.49 


-0.82 


-0.32 


-0.32 


-0.15 


0.01 


-0.27 


-0.12 


-0.3 


-0.2 


0.126 


0.077 


0.058 


0.9 


2.245 


-0.346 


-0.122 


-0.06 


-0.5 


-0.82 


-0.32 


-0.33 


-0.17 


0. 


-0.27 


-0.13 


-0.29 


-0.19 


0.125 


0.077 


0.054 


1.1 


2.435 


-0.349 


-0.115 


-0.05 


-0.52 


-0.82 


-0.33 


-0.35 


-0.19 


-0.02 


-0.27 


-0.15 


-0.28 


-0.19 


0.125 


0.078 


0.05 


1.3 


2.678 


-0.359 


-0.11 


-0.05 


-0.55 


-0.81 


-0.36 


-0.36 


-0.24 


-0.05 


-0.27 


-0.19 


-0.27 


-0.19 


0.127 


0.081 


0.046 


1.5 


2.975 


-0.378 


-0.11 


-0.04 


-0.63 


-0.8 


-0.43 


-0.36 


-0.34 


-0.12 


-0.27 


-0.26 


-0.25 


-0.19 


0.134 


0.086 


0.045 


1.7 


3.329 


-0.409 


-0.119 


-0.04 


-0.85 


-0.75 


-0.78 


-0.32 


-0.65 


-0.31 


-0.31 


-0.44 


-0.22 


-0.18 


0.147 


0.096 


0.053 



TABLE XV: Coefficients for Si{Eq) in the MS scheme as a function of Eq. 



Eo 














5f' 


5f) 


5f) 


^(10) 




^(12) 


^{13) 


^(14) 


^{15) 


^(16) 







1.837 


2.216 


0.729 


0.3 


2. 


-0.31 


1.24 


0.21 


0.43 


-0.26 


1.05 


0.39 


1. 


0.23 


-0.711 


-0.456 


-0.297 


0.5 


1.811 


2.181 


0.715 


0.3 


2.02 


-0.26 


1.26 


0.24 


0.45 


-0.26 


1.06 


0.41 


1. 


0.24 


-0.707 


-0.452 




0.7 


1.775 


2.134 


0.695 


0.29 


2.05 


-0.17 


1.29 


0.3 


0.49 


-0.25 


1.06 


0.43 


1.01 


0.26 


-0.7 


-0.446 




0.9 


1.724 


2.064 


0.664 


0.28 


2.11 


-0.03 


1.34 


0.41 


0.57 


-0.23 


1.08 


0.46 


1.02 


0.29 


-0.691 


-0.437 




1.1 


1.662 


1.971 


0.615 


0.26 


2.21 


0.18 


1.43 


0.59 


0.7 


-0.19 


1.11 


0.53 


1.05 


0.34 


-0.678 


-0.424 




1.3 


1.593 


1.858 


0.542 


0.23 


2.38 


0.5 


1.6 


0.9 


0.92 


-0.12 


1.17 


0.63 


1.09 


0.41 


-0.664 


-0.408 




1.5 


1.532 


1.735 


0.434 


0.16 


2.74 


0.98 


2. 


1.46 


1.38 


0.07 


1.28 


0.85 


1.17 


0.52 


-0.66 


-0.391 




1.7 


1.524 


1.684 0.351 


0.08 


3.76 


1.81 


3.64 


2.88 


2.61 


0.71 


1.59 


1.47 


1.34 


0.72 


-0.698 


-0.396 





once 777.^/2 — £'0 » Aqcd, and are given in the MS scheme by 



ri^^ = 2.1 GeV, 



(2) 



tP = -0.01, 



(8) 



-0.01, 



1 

~~2' 
-0.04, 



.(9) 



-0.02, 
-0.04, 



n 



(5) 



= -0.06, 
0.37, 



and 



^(1) _ rpi2) _ rpi3) _ rp{i) _ ^(6) _ ^(7) _ ^(8) _ ^(13) _ ^(14) _ ^ 

-'2 — -'2 — -'2 — -'2 — -'2 — -'2 — -'2 — -'2 — -'2 — ^1 



(5) 



~12' 



(9) 



-0.04, 



(10) 



-Tl 



(12) 



0.06, 



Tl. 



(11) 



^ = -0.18, 



(A.5) 



(A.6) 



-0.02. 



The remaining, i?o-dependent coefficients of the perturbative corrections are hsted in Ta- 
ble XVII. Since in this case we are expanding the h quark mass about 4.2 GeV, we are 
only showing results for Eq < 2 GeV. The large size of the perturbative corrections to Ti 
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TABLE XVI: Coefficients for 5*2 (-Bq) in the MS scheme as a function of Eq. 



Eo 


5^1) 




si'' 


si'' 


c(5) 


c(6) 
•^2 


0(7) 

•^2 


•^2 


si'' 


si'°' 


0(11) 

•^2 


r.{12) 
•^2 


r.(13) 
•^2 


si''' 


si'" 


o(l(>) 
■^2 


£-,(17) 





0.549 


1.175 


0.78 


0.09 


-5.13 


-1.86 


-1.75 


-3.01 


-6.91 


0.7 


-1.34 


-3.5 


-0.3 


-0.39 


0.085 


-1.169 


-0.187 


0.5 


0.542 


1.16 


0.769 


0.09 


-5. 


-1.8 


-1.65 


-2.91 


-7. 


0.55 


-1.31 


-3.41 


-0.3 


-0.38 


0.164 


-1.131 




0.7 


0.54 


1.155 


0.766 


0.09 


-4.84 


-1.74 


-1.51 


-2.78 


-7.15 


0.35 


-1.27 


-3.3 


-0.29 


-0.36 


0.282 


-1.091 




0.9 


0.544 


1.163 0.774 


0.1 


-4.6 


-1.66 


-1.29 


-2.62 


-7.43 


0.04 


-1.21 


-3.13 


-0.28 


-0.34 0.477 


-1.032 




1.1 


0.554 


1.186 


0.796 


0.12 


-4.28 


-1.57 


-0.97 


-2.43 


-7.9 


-0.37 


-1.14 


-2.9 


-0.28 


-0.32 


0.79 


-0.966 




1.3 


0.567 


1.218 


0.831 


0.15 


-3.88 


-1.46 


-0.48 


-2.19 


-8.7 


-0.9 


-1.05 


-2.61 


-0.27 


-0.29 


1.318 


-0.891 




1.5 


0.57 


1.236 


0.867 


0.19 


-3.44 


-1.34 


0.23 


-1.95 


-10.24 


-1.61 


-0.95 


-2.28 


-0.28 


-0.26 


2.363 


-0.799 




1.7 


0.529 


1.138 


0.786 


0.17 


-3.22 


-1.26 


0.78 


-2. 


-14.43 


-2.79 


-0.94 


-2.07 


-0.33 


-0.25 


5.272 


-0.667 





TABLE XVIL Perturbative coefficients for Ti and T2 in the MS scheme as a function of Eq. 



Eo 


7.(15) 
^1 


7.(16) 
^1 


7.(17) 


7.(15) 
-'2 


7.(16) 
^2 


7.(17) 


1.7 


0.143 


0.083 


-0.009 


0.008 


0.006 


-0.014 


1.8 


0.151 


0.0888 


-0.002 


0.005 


0.004 


-0.013 


1.9 


0.161 


0.095 


0.008 


0.003 


0.003 


-0.011 


2 


0.173 


0.106 


0.03 


0.001 


0.001 


-0.008 



TABLE XVIIL Coefficients for Rq{0,Eo) in the pole scheme as a function of Eq. 





,,'i) 
'\i 






,,'!) 




,,((0 




'1(1 




'\l 


,,(11) 


p'l-:^) 
'!() 




,,(14) 




,,("0 


,,'17) 

''.) 


0.5 


0.973 


-0.002 


-0.001 


0. 


0. 


-0.01 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0.7 


0.93 


-0.004 


-0.002 


0. 


-0.01 


-0.02 


0. 


-0.01 


0. 


0. 


0. 


0. 


0. 


0. 


0.001 


0. 


0. 


0.9 


0.86 


-0.009 


-0.004 


0. 


-0.01 


-0.04 


-0.01 


-0.02 


0. 


0.01 


-0.01 


0. 


-0.01 


-0.01 


0.001 


0. 


0. 


1.1 


0.761 


-0.016 


-0.007 


0. 


-0.02 


-0.08 


-0.02 


-0.03 


-0.01 


0.01 


-0.01 


0. 


-0.01 


-0.01 


0.001 


0. 


0. 


1.3 


0.634 


-0.025 


-0.01 


0. 


-0.03 


-0.11 


-0.02 


-0.05 


-0.01 


0.02 


-0.01 


0. 


-0.02 


-0.02 


0.001 


-0.001 


0. 


1.5 


0.483 


-0.038 


-0.014 


-0.01 


-0.05 


-0.15 


-0.03 


-0.07 


-0.01 


0.02 


-0.02 


0. 


-0.03 


-0.03 


0. 


-0.002 


-0.001 


1.7 


0.318 


-0.054 


-0.018 


-0.01 


-0.06 


-0.19 


-0.04 


-0.08 


-0.01 


0.02 


-0.02 


-0.01 


-0.03 


-0.03 


-0.001 


-0.003 


-0.002 



(compared to its values in the l^* or PS schemes) occur to try to compensate for the bad 
choice of mass scheme. 



4. The pole mass scheme 



The B — > Xciu decay width in the pole scheme is given by 

\V l"^ r 

r{B XJP) = m^s 0.370 - 0.115 A - 0.012 A^ + 0. A=^ 

19277"^ 

-0.04 Ai - 0.10 A2 - 0.01 AiA + 0.02 A2A 
-0.02 pi + 0.02 p2 - 0.02 Ti + 0. T2 

-0.03 Tg- 0.02 7;- 0.040 e- 0.022 e^LM + 0.007 eA , (A.7) 

We tabulate the shape variables defined in Eq. (1) in Tables XVIII, XIX, and XX, and 
those defined in Eq. (2) in Tables XXI and XXII in the pole mass scheme. For Si and S2 
we do not show the £'o-dependence of the order eA terms, as they are not known. For all 
quantities the coefficients of the expansions are defined as in Eq. (15). 

For the B — > shape variables defined in Eq. (4), tI;^^\ . . . T^^^^^ are independent of Eq, 



20 



TABLE XIX: Coefficients for Ri(Eq) in the pole scheme as a function of Eq. 



Eq 


rY' 


r[' 


R^' 


RY^' 


Rf> 


Rf> 


rY' 


RY' 


RY' 


RY ' 


RY ' 


("121 

RY ' 


RY ' 


C141 

RY ' 


RY ' 


RY ' 


RY ' 





1.429 


-0.054 


-0.014 


0. 


-0.07 


-0.18 


-0.04 


-0.05 


-0.03 


0.01 


-0.02 


-0.01 


-0.03 


-0.03 


0. 


-0.002 


-0.001 


0.5 


1.459 


-0.054 


-0.014 


0. 


-0.07 


-0.18 


-0.03 


-0.05 


-0.03 


0.01 


-0.02 


-0.01 


-0.03 


-0.03 


0. 


-0.002 


-0.001 


0.7 


1.498 


-0.054 


-0.013 


0. 


-0.07 


-0.17 


-0.03 


-0.05 


-0.03 


0.01 


-0.02 


-0.02 


-0.03 


-0.03 


-0.001 


-0.003 


-0.001 


0.9 


1.554 


-0.054 


-0.013 


0. 


-0.07 


-0.17 


-0.03 


-0.05 


-0.04 


0.01 


-0.02 


-0.02 


-0.03 


-0.03 


-0.001 


-0.003 


-0.001 


1.1 


1.625 


-0.055 


-0.012 


0. 


-0.07 


-0.16 


-0.03 


-0.05 


-0.04 


0. 


-0.02 


-0.02 


-0.02 


-0.03 


-0.002 


-0.003 


-0.001 


1.3 


1.71 


-0.056 


-0.012 


0. 


-0.08 


-0.15 


-0.03 


-0.05 


-0.04 


0. 


-0.02 


-0.02 


-0.02 


-0.03 


-0.002 


-0.003 


-0.001 


1.5 


1.806 


-0.058 


-0.013 


0. 


-0.08 


-0.14 


-0.04 


-0.05 


-0.05 


-0.01 


-0.02 


-0.03 


-0.02 


-0.03 


-0.002 


-0.003 


-0.001 


1.7 


1.913 


-0.06 


-0.013 


0. 


-0.1 


-0.13 


-0.05 


-0.04 


-0.07 


-0.02 


-0.02 


-0.04 


-0.02 


-0.02 


-0.003 


-0.003 


-0.001 



TABLE XX: Coefficients for R2{Eq) in the pole scheme as a function of Eq. 





^2 




^2 










^2 


r(9) 
^2 


^2 




^2 


"■2 






^2 







2.241 


-0.184 


-0.041 


-0.01 


-0.26 


-0.58 


-0.11 


-0.15 


-0.13 


0.02 


-0.08 


-0.06 


-0.09 


-0.11 


-0.002 


-0.008 


-0.003 


0.5 


2.299 


-0.185 


-0.041 


-0.01 


-0.26 


-0.58 


-0.11 


-0.15 


-0.14 


0.02 


-0.08 


-0.06 


-0.09 


-0.11 


-0.003 


-0.009 


-0.003 


0.7 


2.388 


-0.188 


-0.04 


-0.01 


-0.26 


-0.57 


-0.11 


-0.15 


-0.14 


0.02 


-0.08 


-0.07 


-0.09 


-0.11 


-0.004 


-0.01 


-0.003 


0.9 


2.529 


-0.193 


-0.04 


-0.01 


-0.26 


-0.57 


-0.11 


-0.15 


-0.15 


0.01 


-0.08 


-0.07 


-0.09 


-0.11 


-0.005 


-0.011 


-0.003 


1.1 


2.726 


-0.2 


-0.04 


-0.01 


-0.27 


-0.56 


-0.11 


-0.15 


-0.16 


0. 


-0.08 


-0.08 


-0.09 


-0.11 


-0.007 


-0.012 


-0.003 


1.3 


2.981 


-0.211 


-0.041 


-0.01 


-0.29 


-0.55 


-0.12 


-0.15 


-0.18 


-0.02 


-0.08 


-0.1 


-0.08 


-0.1 


-0.008 


-0.012 


-0.003 


1.5 


3.298 


-0.225 


-0.043 


-0.01 


-0.33 


-0.54 


-0.14 


-0.15 


-0.22 


-0.04 


-0.09 


-0.12 


-0.08 


-0.1 


-0.01 


-0.013 


-0.003 


1.7 


3.678 


-0.243 


-0.047 


-0.01 


-0.4 


-0.52 


-0.19 


-0.14 


-0.32 


-0.09 


-0.1 


-0.17 


-0.07 


-0.1 


-0.01 


-0.013 


-0.003 



TABLE XXI: Coefficients for Si{Eo) in the pole scheme as a function of Eq. 



Eo 




5f) 






5f) 


5f) 


5f) 


5f) 




^(10) 




^(12) 


^(13) 




^(15) 


^(16) 










1.248 


0.262 


0.06 


1.02 


-0.32 


0.41 


-0.11 


0.42 


-0.21 


0.3 


0.15 


0.28 


0.08 


0.102 


0.111 


0.038 


0.5 





1.231 


0.258 


0.06 


1.03 


-0.29 


0.41 


-0.09 


0.43 


-0.21 


0.31 


0.16 


0.28 


0.08 


0.097 0.107 




0.7 





1.209 


0.251 


0.06 


1.05 


-0.25 


0.42 


-0.07 


0.45 


-0.21 


0.31 


0.17 


0.28 


0.09 


0.092 


0.102 




0.9 





1.18 


0.241 


0.06 


1.08 


-0.18 


0.44 


-0.03 


0.48 


-0.2 


0.31 


0.18 


0.29 


0.1 


0.084 


0.095 




1.1 





1.148 


0.228 


0.05 


1.13 


-0.08 


0.46 


0.03 


0.54 


-0.19 


0.32 


0.21 


0.29 


0.12 


0.075 


0.086 




1.3 





1.118 


0.211 


0.04 


1.22 


0.04 


0.51 


0.12 


0.63 


-0.16 


0.34 


0.26 


0.3 


0.14 


0.064 


0.077 




1.5 





1.1 


0.194 


0.04 


1.38 


0.2 


0.6 


0.26 


0.8 


-0.11 


0.37 


0.35 


0.31 


0.17 


0.054 


0.067 




1.7 





1.112 


0.188 


0.03 


1.7 


0.4 


0.83 


0.47 


1.16 


0.04 


0.43 


0.53 


0.32 


0.21 


0.044 


0.057 





once uib/^ — £'o ^ Aqcd, and are given in the pole mass scheme by 

_ t(2) _ 1 ^(3) _ ^(4) _ ^(5) _ ^(7) _ 

^1 ~ ~2~' ~ 2' "-'i "-'i — -'-1 — 

rf ^ = -0.14, Ti^'^ = -0.03, Ti^^ = -0.02, T^^^^ = 0.11, 

yUi) ^ r^(m ^ 0., Ti^""^ = Ti^^^ = -0.03, (A.8) 

and 

rpW ^ 2^(2) ^ rp{3) ^ rp(4.) ^ rp(6) ^ rp(7) ^ rp(8) ^ ^(13) ^ ^(14) ^ q ^j^ g. 

T^^^ = ^ = -0.03, Ti^°^ = -T^^^^ = 0.05, T^^^^ = -0.02. 

The remaining, E'o-dependent coefficients of the perturbative corrections are hsted in Ta- 
ble XXIII. 
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TABLE XXII: Coefficients for 5*2 (-Eg) in tlie pole scheme as a function of Eq. 



Eo 




sr 




sr 


sr 


sr 




sr 


sr 


sr' 


sr' 


sr' 


sr' 


sr' 


sr' 


sr' 


sr' 











0.297 


0.1 


-3.9 





-0.86 


-0.82 


-4.52 


1.24 


-0.73 


-2.2 








0.301 


0.255 


0.146 


0.5 








0.294 


0.1 


-3.84 





-0.83 


-0.8 


-4.57 


1.16 


-0.72 


-2.17 








0.273 


0.235 




0.7 








0.293 


0.1 


-3.77 





-0.78 


-0.79 


-4.66 


1.04 


-0.71 


-2.13 








0.241 


0.212 




0.9 








0.296 


0.1 


-3.67 





-0.71 


-0.77 


-4.83 


0.87 


-0.69 


-2.07 








0.202 


0.182 




1.1 








0.301 


0.1 


-3.56 





-0.61 


-0.75 


-5.1 


0.65 


-0.67 


-2.01 








0.16 


0.149 




1.3 








0.307 


0.11 


-3.46 





-0.48 


-0.72 


-5.56 


0.37 


-0.65 


-1.95 








0.12 


0.115 




1.5 








0.306 


0.11 


-3.4 





-0.34 


-0.69 


-6.39 


0.02 


-0.64 


-1.92 








0.083 


0.084 




1.7 








0.287 


0.1 


-3.43 





-0.27 


-0.62 


-8.05 


-0.43 


-0.65 


-1.94 








0.051 


0.056 





TABLE XXIII: Perturbative coefficients for Ti and T2 in the pole scheme as a function of Eq. 



Eo 


^^(15) 


^1 


^(17) 


^(15) 




71(17) 
-'2 


1.7 


-0.077 


-0.069 


0.008 


0.022 


0.014 


-0.007 


1.8 


-0.074 


-0.068 


0.012 


0.02 


0.013 


-0.009 


1.9 


-0.071 


-0.067 


0.016 


0.017 


0.012 


-0.011 


2 


-0.068 


-0.065 


0.021 


0.015 


0.011 


-0.012 


2.1 


-0.063 


-0.062 


0.026 


0.012 


0.009 


-0.013 


2.2 


-0.058 


-0.059 


0.031 


0.009 


0.008 


-0.013 
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